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Approximating electronically excited states with equation-of-motion linear coupled-cluster theory 
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A new perturbative approach to canonical equation-of-motion coupled-cluster theory is presented using 
coupled-cluster perturbation theory. A second-order Mpller-Plesset partitioning of the Hamiltonian is used 
to obtain the well known equation-of-motion many-body perturbation theory (EOM-MBPT(2)) equations 
and two new equation-of-motion methods based on the linear coupled-cluster doubles (EOM-LCCD) and 
linear coupled-cluster singles and doubles (EOM-LCCSD) wavefunctions. This is achieved by performing 
a short-circuiting procedure on the MBPT(2) similarity transformed Hamiltonian. These new methods are 
benchmarked against very accurate theoretical and experimental spectra from 25 small organic molecules. It 
is found that the proposed methods have excellent agreement with canonical EOM-CCSD state for state or¬ 
derings and relative excited state energies as well as acceptable quantitative agreement for absolute excitation 
energies compared with the best estimate theory and experimental spectra. 


I. INTRODUCTION 

Optical spectroscopy is an important and ubiquitous 
tool in modern experimental studies. Improvements over 
the years in the use of lasers in optical spectroscopy has 
brought the field to an impressive level of accuracy and 
precision that is difficult to match using modern theoret¬ 
ical methods. Even so, theory has an important role to 
play in optical spectroscopy. By assisting in excited state 
assignments and predicting the density of excited states 
expected in a given energy range, computational results 
can be a powerful tool for spectroscopists. 

Widespread use of time dependent density functional 
theory^ (TD-DET) for the calculation of excited states 
will continue for a long time to come due to the attrac¬ 
tive low computational scaling inherent to the method. 
However, reliability and accuracy problems^ in TD-DFT 
keep the theory from replacing the more computation¬ 
ally expensive but highly accurate multireference con¬ 
figuration interaction^ (MRCI) or equation of motion 
coupled-clusteri^ (EOM-CC) theories. The most com¬ 
monly used form of EOM-CC theory is the EOM-CCSD 
variant)^ which limits ground and excited state excita¬ 
tions to singles and doubles only. This approximation is 
qualitatively consistent and approaches quantitative ac¬ 
curacy in many cases where single excitations are dom¬ 
inant. For doubly excited states, perturbative^ (EOM- 
CCSD(T)) or iterative^ (EOM-CCSDT-n) triple excita¬ 
tions are required^ii^ to obtain quantitative agreement 
with experiment. 

The development of approximate excited state meth¬ 
ods based on many-body perturbation and coupled- 
cluster theory has been an active and fruitful ende- 
vour since the emergence of the heldiii Many ap¬ 
proaches make use of a mean-field starting reference 
based on configuration interaction singles^^ and add per- 
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turbative corrections to that referenceJ^^— Other ap¬ 
proaches use a many-body treatment of propagator the¬ 
ory such as the algebraic-diagrammatic constructio n^^i^° 
(ADC) scheme or by using approximate coupled-cluster 
linear-response2U— theory with the CCn methods<^ Ap¬ 
proximate methods based on equation-of-motion the¬ 
ory can either add a post hoc EOM-CC perturbative 
correctio nb^^i^^ for triples or an a priori perturbative 
approximation.—^— We will use perturbation theory here 
to develop a new perturbative EOM-CC method that also 
includes infinite-order effects in the ground state wave- 
function. The goal will be to provide a consistent and 
accurate way to obtain excitation energies from either 
a linear CCD (LCCD) or CCSD (LCCSD) ground state 
wavefunction. These methods are not intended to be a re¬ 
placement method for existing fast (EOM-MBPT(2)) or 
accurate (EOM-CCSD) methods when considering just 
the calculation of excited states, but to provide a route 
toward a consistent excited state spectra from a specific 
ground state wavefunction. 

In this work, we will use coupled-cluster perturbation 
theory^S (CCPT) to derive a general EOM-CCPT frame¬ 
work in Sec. im The various EOM-CC approximations 
will then be systematically derived from EOM-CCPT 
by a choice in the Hamiltonian perturbation partition¬ 
ing. Numerical results of these approximations compared 
against very accurate theoretical and experimental spec¬ 
tra are presented in Sec. El 


II. THEORY 

A. Coupled-Cluster Theory 

Central to coupled-cluster theory^ is the similarity 
transformation of the electronic Hamiltonian by the ex¬ 
ponential wave operator 

H = e-^He^. (1) 
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The excitation cluster operators (here limited to single 
and double excitations) are given by 

f = fi + ^ ttan + i ^ (2) 

ai abij 

and the normal ordered electronic Hamiltonian is 


Here uik is the fc’th general eigenvalue, and R{k) {L{k)) 
is the fc’th right (left) general eigenvector solution to the 
non-Hermitian matrix, H. The {)c denotes that only 
fully connected contractions are included. These eigen¬ 
vector solutions are defined as the linear excitation and 
de-excitation operators (limited again to only singles and 
doubles) 


H = {(j>o\H\(j)o) + f^{pU} + < 1 ^ sr} (3) 

pq pqrs 

= {cl)o\H\cj)o)+F + W. (4) 


Here / is the one-electron Fock matrix, v the antisym¬ 
metrized two-electron integrals, } denotes normal or¬ 
dering of the enclosed operators and the one- {F) and 
two-particle (W) are defined appropriately. Throughout 
we reserve the indices i,j, k, - ■ ■ and a,b,c, - ■ ■ for occu¬ 
pied and virtual orbitals respectively while the indices 
p^q, - ■ ■ may refer to either occupied or virtual orbitals. 
Because the Hamiltonian only contains one- and two- 
particle operators, the Baker-Campbell-Hausdorff expan¬ 
sion of Eq. [T] naturally truncates after four commutators 
giving 


H = H 

1 rrr. 
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Starting with an appropriate single-reference mean- 
field reference, |(()o), the coupled-cluster Schrbdinger 
equation is expressed with Eq. [5] as 

H\(l,o) = Ecc\M- ( 6 ) 

The CC correlation energy given by projecting on the left 
by the reference function 


R{k) = ro(fc) -I- Ri{k) + R 2 (k) = 

ai abij 


L{k) = Li{k)+L2{k) =Y^lik)i^^+^Yiaf^{k)i'<apb, 

ai abij 

( 12 ) 

with the biorthogonalization constraint 

{ML{k)R{k'Mo) = hk'- (13) 

Note that because R{k) is an excitation operator, it will 
commute with the ground state cluster amplitudes 


f,i?(fc) 


= 0 . 


(14) 


B. Coupled-Cluster Perturbation Theory 
Similarity Transformed Hamiltonian 

As is always possible in Hamiltonian based perturba¬ 
tion theory we can represent the complete Hamiltonian 
in terms of a zero’th order Hamiltonian F[q and pertur¬ 
bation V 


H = Ho+ aV 


(15) 


{cl)o\H\M = Fee (7) 

while the necessary cluster amplitudes are obtained by 

solving the equations generated by projecting on the left 
with the auxiliary space 

(</),|i7|4)=0 (8) 

where {(pg \ is the set of 5 -fold excited determinants 

(<(>g| = (<(>ol4hi4«2 ■ •-* 3 % (9) 


where a —>■ 1 is an order parameter used for convenience. 
Equation [ 15 ] can be directly inserted into in any coupled- 
cluster similarity transformation or energy functional and 
expanded to the desired order in a. Examples other 
than the usual form given by Eq. [5] would include the 
Hermitia n^^i^^ (e^^f/e^) or unitary^ {e^"' He^) expan¬ 
sions, but they do not simply terminate without trunca¬ 
tion. 

To facilitate the decomposition into a reference and 
perturbation operator, we rearrange the Hamiltonian 
into particle excitation rank form 


The conventional approach to computing excited states 
through the equation-of-motion (EOM) method;^ which 
directly computes the fc’th excitation energy (wfc) relative 
to the ground state coupled-cluster wavefunction, is also 
obtained with Eq. [5] as 


{(j)o\L{k) H,R{k) \(l)o)c = uJk. 


( 10 ) 


H= (16) 

where the [n] superscript denotes that the operator 
changes particle rank from right to left by n. Erom this 
point the CC perturbation framework used is the same 
as we have employed in the pasti ^^’^°’^^ Here the expo¬ 
nential wavefunction is expanded to the appropriate 
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order in terms of n’th order cluster operators The 

computation of these perturbative cluster operators is 
done order by order using standard perturbation theory. 

I 


This is easily illustrated by expanding Eq. [5] using the 
partitioned Hamiltonian of Eq. 1151 The similarity trans¬ 
formed Hamiltonian to arbitrary order in a is given by 


[(^HoSn,n' + + ^ ^ + VSr,,n'+m+l) 
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? E [[[{HoSn,n'+m+m' +V6r,,n'+m+m' + l) , f <”') 




n'mm'm" 


Throughout this work we denote the exclusive order in a 
of an operator with the superscript (n) and we use {n} 
to denote an expansion inclusive of all orders in a up to 
n. With this notation, the general n’th order similarity 
transformed Hamiltonian can be expressed as a compact 
series 


ijl"} = iJp-h • • •-b (18) 

The n’th order ground state CCPT cluster amplitudes 
are completely defined from insertion of the exclusive 
similarity transformed Hamiltonian, Eq. 1171 into [8] 

(</),|i7('^)|</)o)=0. (19) 

Using these n’th order ground state amplitudes in Eq. 
[7] will give the n -|- 1 order correlation energy beginning 
with n = 2 


( 20 ) 

From this the exact coupled-cluster correlation energy, 
relative to the mean-field energy 

AEcc=Ecc-{0\H\0), (21) 

is 

AEcc - AE^^I + AE^^l + • • • + ( 22 ) 

or equivalently 

AEcc - + • • • + (23) 

= (24) 


To obtain excitation energies relative to the n’th order 
energy without perturbatively expanding the L and R op¬ 
erators, we use the n’th order inclusive similarity trans¬ 
formed Hamiltonian (Eq. ITS)) in Eq. [10] rather than the 
exclusive Eq. [171 giving 


{ML{k) fi7^”>,i?(fc)l l^o) =c^fe- 


(25) 


There are many ways to partition the Hamiltonain. We 
have used the particle rank conserving partitioning of the 
Hamiltonian, 


Hq = (0|77|0) -b f[o] + bUM (26) 

V = -b , (27) 

in previous wor k^^i^°i^^ to great success. Alternatively 
we can also choose the traditional Mpller-Plesset (MP) 
partitioning which we will discuss next. 


C. Second-Order M0ller-Plesset Partitioning 


The standard approach in electronic structure the¬ 
ory is to use generalized many-body perturbation the¬ 
ory (GMBPT)^ based on the Mpller-Plesset (MP) par¬ 
titioning of the Hamiltonian, given by 

Ho = (28) 

V = E^^^i + WM -b -b . (29) 

Using Eqs. (281 and (291 in Eq. [T3 the first two orders of 
Hin) 

are respectively 


^(1) 




+ U 


(30) 


and 


ij(2) 






(31) 


Inserting Eq. |30| into Eq. [T9| gives the quintessential 
GMBPT(l) wavefunction 

(<(.i|F[olfi^^)+F[+il|((.o)c = 0 (32) 

+ W[+2l|((,o)c = 0 (33) 

from which 

U(2) = + W^-^^f^'^^\cl>o)c- (34) 
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Ethylene* Butadiene* Hexatriene Octatetraene 



Benzene Naphthalene Furan Pyrrole Imidazole 



Formaldehyde* Acetaldehyde* Acetone* Benzoquinone 



Formamide Acetamide Propanamide 


FIG. 1. Member molecules of the Miilheim test set. Names 
denote with an * are also members of the Yale set, with Ac¬ 
etaldehyde only a member of the Yale set. 


The GMBPT(2) wavefunction is 


^ -b + (#[-1] -b Vo)c = 0 

(35) 

= 0. (36) 

The triples contribution does not contribute to the first- 
or second-order wavefunction. The inclusive second- 
order similarity transformed Hamiltonian, 


fjm 


-bH-b 


^[0] ^ y(l) 


plo] ^ ^(2) 




(37) 

is the GMBPT generalization of previous^! ^ MBPT ap¬ 
proximations to H. The spin-orbital equations for the 
matrix elements^li^ using Einstein notation are 


n] =f; - fitf^ - , (38) 

ni =K + , (39) 

(40) 

n: =0, (41) 


-yab —f.ab , j-amMl) , }_;p,mn,ab(l) 

^cd ^cd ' ^cd ' 2^^^ ^'fnn ■> 

(42) 

'^kt 

(43) 

=< - + P{^J)P{ab)v2^ 

,ae{l) 
^mj ’ 

(44) 

'^ka =^L-^e'atk^\ 

(45) 

'L/ai 1 ;r,mi^a{l) 

^bc ^bc ' ^bc bjn i 

njia _;ria , -:im ,a(l) -^io,e(l) . ni,ea{l) 

^jk ^jk ' ^jk ^je^k ' Je^jk ' 

(46) 

2^e/Gfc bmk > 

(47) 

njab —r.ab _ -af).e(l) , -amMl) _ 

‘^ci ^ci ^ce'^i ' ^ci 



(48) 

Kb =^b, 

(49) 

K =0, 

(50) 

■TjO-jb _-mj,ab{l) 

^cdi ^cd ^mi ’ 

(51) 

'T-fcbf ~ ’ 

(52) 

Kfk =7^(«&)Ci?'^ 

(53) 

where P{pq) is the anti-permutation operator 


P{pq)9{pq ■■■rs) = g{pq ...rs)- g{pq ...sr). 

(54) 


The difference between GMBPT and MP theory occurs 
when /“, fj and f§ p 0 as they would be for canoni¬ 
cal HE orbitals. The first introduces a first-order singles 
wavefunction (Eq. 1321) , while the second and third terms 
have to be summed to all order as in coupled-cluster 
theory or better transformed away by a semi-canonical 
occupied-occupied and virtual-virtual rotation. The first- 
order t amplitudes used in Eqs. [38ll5^ are obtained from 
Eq. [22] and Eq. [321 Second-order amplitude contribu¬ 
tions are limited in this case to Hf and because in 
the GMBPT partitioning can only give particle 

excitation matrix elements. Because these matrix ele¬ 
ments are by definition zero (Eq. jSj) there are no actual 
contributions from the second-order amplitudes to the 
excited state spectrum. 


D. Inclusive Second-Order M0ller-Plesset 
Partitioning 

The EOM-MBPT(2) approach theoretically suffers 
from the fact that the starting wavefunction is strictly 
a single shot perturbation calculation, while the diago- 
nalization of H introduces infinite-order contributions to 
the excitation wavefunction. This can be overcome, while 
using the GMBPT partitioning, by short-circuiting the 
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FIG. 2. Comparison against the Miilheim data set EOM-CCSD values. 


order-by-order amplitude computation to instead solve 

(</>5|ff{2}|0o)=o. (55) 

Here we combine the first- and second-order amplitude 
equations to obtain 

^[-i]^aoo) ^[+1] 1^^^^ ^ Q 

EauationsISSlandlSTlare incidentally the linearized CCSD 
amplitude equations (LCCSD). If the singles contribu¬ 
tions in Eq. [56]are removed, the resulting reference wave- 
function is then simply linearized CCD (LCCD).— This is 
a convenient approximation as it requires no new terms in 
the similarity transformed Hamiltonian, while introduc¬ 
ing infinite-order character to the reference wavefunction. 


III. ELECTRONIC STRUCTURE 
CALCULATIONS 

To assess both the quality and general behavior of 
these EOM approximations when computing vertical ex¬ 
citation energies, we employ two established sets of gas 
phase molecules (see Figure [T]): the set of 24 organic 
molecules from Schreiber et (referred to hereafter 

as the Miilheim set) and the set of 11 organic molecules 
from Caricato et a/.— (referred to here as the Yale set). 
This Miilheim set contains 121 reference single excitation 
energies^! which use molecular geometries obtained with 
an MP2/6-31G* optimization in which the first row Is 
atomic orbitals are dropped. The extent and balanced 
nature of the Miilheim set readily lends itself to obtain¬ 
ing statistics between different theoretical methods. To 


State MBPT(2) 

LCCD 

LCCSD 

CCSD 

Exp. 

ethylene 

B^u 

7.15 

7.43 

7.45 

7.28 

7.11 

Biu 

7.76 

8.13 

8.18 

7.97 

7.65 

Big 

7.79 

8.08 

8.10 

7.97 

7.80 

B2g 

7.84 

8.12 

8.14 

7.93 

7.90 

Ag 

8.19 

8.46 

8.48 

8.31 

8.28 

Bsu 

8.64 

8.92 

8.94 

8.77 

8.62 

Bsu 

8.93 

9.20 

9.22 

9.05 

8.90 

Bsu 

9.05 

9.33 

9.35 

9.17 

9.08 

Big 

9.19 

9.46 

9.49 

9.34 

9.20 

Biu 

9.18 

9.47 

9.50 

9.32 

9.33 

Bsu 

9.84 

10.11 

10.13 

9.96 

9.51 

isobutene 

Bi 

6.29 

6.54 

6.56 

6.38 

6.17 


6.81 

7.06 

7.10 

6.91 

6.70 

trans-1,3-butadiene 





Bu 

6.12 

6.54 

6.63 

6.29 

5.91 

Bg 

6.17 

6.47 

6.49 

6.24 

6.22 

Au 

6.47 

6.79 

6.80 

6.55 

6.66 

Bu 

7.06 

7.39 

7.43 

7.16 

7.07 

Bg 

7.28 

7.59 

7.60 

7.36 

7.36 

Ag 

7.54 

7.85 

7.85 

7.61 

7.62 

Bu 

8.08 

8.39 

8.40 

8.15 

8.00 

MAD 

0.09 

0.31 

0.34 

0.14 


RMS 

0.12 

0.34 

0.37 

0.19 


MD 

0.33 

0.63 

0.72 

0.45 



TABLE I. Excitation energies for a few small alkene molecules 
computed with the d-aug-cc-pVDZ basis using MP2/6- 
311-l-G** geometries. Units are in electron volts (eV), ex¬ 
perimental references can be found in Ref. Q. 


compare the approximate EOM-CC methods proposed 
here, we take our best theoretical reference to be EOM- 
CCSDT-3)^ which is known to give both accurate re¬ 
sults and is a systematic EOM-CC type theory which will 
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FIG. 3. Comparison against the Miilheim data set EOM-CCSDT-3 values. 


provide a balanced comparison. The Yale set choice in 
molecules is a subset of the Miilheim set with the simple 
addition of acetaldehyde, but with a different selection 
in excitation energies (69 in all) whose reference values 
are taken from gas phase experiments (the experiment 
literature references can be found in Ref. d). This com¬ 
parison to experiment provides a different and valuable 
perspective with which to augment theory vs. theory 
examinations. Geometries for the Yale set are obtained 
from an MP2/6-311-fG** optimization. 


Galculations using the Miilheim set typically^ use the 
TZVP basis set (3slp on hydrogen 5s3pld on first row 
atoms) from Schafer et al^ This is a light weight ba¬ 
sis set that allows the use of many different kinds of 
computational methods with only a modest cost in over¬ 
all accuracyi^ While perfectly acceptable when compar¬ 
ing different methods to each other, this basis is not 
accurate enough to compare with experimental results. 
This is especially true when calculating Rydberg states, 
where the inclusion of many diffuse functions is keyj ^^d^ 
When computing the excitation energies for the Yale 
set, we use the doubly diffuse function augmented Dun¬ 
ning basis )^^d6 d-aug-cc-pVDZ (4s3p on hydrogen and 
5s4p3c? on first row atoms). This basis is comparable to 
the 6-311(3-1-,3-|-)G** basis used by Caricato et al. in 
their calculations^ d^d^ with almost no gain in including 
a third set of diffuse functions4^ 


All of the reported electronic structure results were ob¬ 
tained using the serial AGESII^ and parallel Aces4^ ab 
initio quantum chemistry packages. Calculations were 
performed on the University of Florida HiPerGator high 
performance cluster. Throughout this paper we use RMS 
to mean ’’root mean square,” MAD to mean ’’mean av¬ 
erage absolute deviation,” and MD to mean ’’absolute 
maximum deviation.” 


State 

MBPT(2) 

LCCD 

LCCSD 

CCSD 

Exp. 

acetaldehyde 

A" 4.19 

4.28 

4.57 

4.32 

4.28 

A' 

6.64 

6.68 

6.94 

6.78 

6.82 

A' 

7.32 

7.59 

7.61 

7.67 

7.46 

A' 

7.56 

7.37 

7.82 

7.46 

7.75 

A' 

8.24 

8.27 

8.51 

8.36 

8.43 

A' 

8.27 

8.31 

8.54 

8.39 

8.69 

acetone 

A 2 

4.38 

4.44 

4.75 

4.48 

4.43 

B 2 

6.25 

6.29 

6.54 

6.37 

6.36 

A 2 

7.17 

7.19 

7.44 

7.27 

7.36 

Al 

7.26 

7.29 

7.55 

7.37 

7.41 

B 2 

7.25 

7.27 

7.51 

7.35 

7.49 

Al 

7.88 

7.90 

8.15 

7.98 

7.80 

B 2 

7.65 

7.68 

7.93 

7.76 

8.09 

Bi 

7.95 

7.98 

8.22 

8.06 

8.17 

formaldehyde 

A 2 3.82 

3.96 

4.23 

3.99 

4.00 

B 2 

6.87 

6.93 

7.17 

7.03 

7.08 

B 2 

7.70 

7.75 

7.96 

7.83 

7.97 

Al 

7.83 

7.87 

8.10 

7.97 

8.14 

A 2 

8.06 

8.09 

8.32 

8.19 

8.37 

B 2 

8.77 

8.82 

9.05 

8.91 

8.88 

Bi 

9.10 

9.14 

9.05 

9.23 

9.00 

A 2 

9.20 

9.24 

9.45 

9.32 

9.22 

B 2 

9.03 

9.07 

9.29 

9.16 

9.26 

Al 

9.07 

9.12 

9.35 

9.20 

9.58 

B 2 

9.07 

9.10 

9.32 

9.19 

9.63 

MAD 

0.22 

0.18 

0.14 

0.13 


RMS 

0.26 

0.23 

0.17 

0.18 


MD 

0.56 

0.53 

0.35 

0.44 


TABLE 

II. Excitation energies 

! for a few 

small carbonyl 

molecules computed 

with the 

d-aug-cc-pVDZ basis 

using 


MP2/6-311-I-G** geometries. Units are in electron volts (eV), 
experimental references can be found in Ref. 0. 
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IV. NUMERICAL RESULTS AND 
DISCUSSION 

A. EOM-MBPT(2) 

Developed two decades ago for canonical HF— EOM- 
MBPT(2) is a well established perturbation method 
which has been largely ignored by the community in fa¬ 
vor of other approximations such as CC32i or ADC(2)Fi^ 
The EOM-MBPT(2) method remains relevant because 
of the non-iterative O(n^) computational cost of the 
ground-state wavefunction and the general accuracy of 
MBPT(2). It is also reasonably accurate for equilibrium 
vertical excitation energies. This accuracy, within the 
Miilheim test set, is illustrated in the correlation fig¬ 
ures [^a), EJa) and Table HVl Compared to the standard 
EOM-CCSD method, EOM-MBPT(2) has an average de¬ 
viation of 0.13 eV while comparison with EOM-CCSDT-3 
gives an average deviation of 0.31 eV respectively, con¬ 
sistent with the known ^ 0.2 eV triples correction to 
EOM-CCDSD. 

The EOM-MBPT(2) approximation and its Lowdin 
partitioned varian t^^’^^ was recently benchmarked (us¬ 
ing the 6-311(3-1-,3-l-)G** basis set) against EOM-CCSD 
and experiment using the Yale test set by Goings et al^ 
Our computed EOM-MBPT(2)/d-aug-cc-pVDZ values, 
presented in Tables III m and uni yield nearly identi¬ 
cal results with minor differences due to the alternative 
basis set choice. The EOM-MBPT(2) method is com¬ 
pletely satisfactory when compared to experimental Ry¬ 
dberg states (such as with the alkenes and carbonyl’s), 
but is much less reliable for the n —> tt* and tt —>■ tt* 
heterocycle valence states. The comparative difficulty in 
describing valence states is a limitation of EOM theory in 
general which is only satisfactorily overcome in CC the¬ 
ory by using methods that include triples or other meth¬ 
ods such as similarity transformed equation-of-motion 
coupled-cluster theory^-^ (STEOM-CC). 


B. Linear EOM-CC 

By employing Eqs. [50] and |57| to compute the am¬ 
plitudes used in the second-order similarity transformed 
Hamiltonian (Eq. l37l) . we now have an equation-of- 
motion theory completely consistent with with the LCCD 
(LCCSD) wavefunction that only includes terms in the 
similarity transformed Hamiltonian up to while en¬ 
suring 

{^g\H\^o) = 0. (58) 

These linear EOM methods have the same iterative 
0{n^) computational scaling as canonical EOM-CCSD. 
However, with the removal of the quadratic amplitude 
intermediates, the computational overhead and I/O de¬ 
mands are significantly reduced. Benchmark data from 
the Miilheim set shows a systematic overestimate of exci- 


State MBPT(2) 

LCCD 

LCCSD 

CCSD 

Exp. 

pyrazine 

Bsu 

4.52 

4.63 

4.63 

4.33 

3.83 

B2u 

5.36 

5.51 

5.50 

5.10 

4.81 

B2g 

6.17 

6.26 

6.29 

6.01 

5.46 

Big 

7.14 

7.35 

7.34 

7.07 

6.10 

Biu 

7.04 

7.28 

7.33 

6.95 

6.51 

pyridazine 

Bi 

4.03 

4.32 

4.36 

4.03 

3.60 


5.33 

5.33 

5.75 

5.33 

5.00 

^2 

4.63 

4.63 

4.96 

4.63 

5.30 

Bi 

6.62 

6.62 

6.93 

6.62 

6.00 

B2 

6.26 

6.26 

6.43 

6.26 

6.50 

pyridine 

Bi 

5.28 

5.38 

5.46 

5.17 

4.59 

B2 

5.41 

5.57 

5.61 

5.23 

4.99 

2^2 

5.64 

5.80 

5.92 

5.60 

5.43 


6.76 

6.86 

6.92 

6.69 

6.38 

pyrimidine 

Bi 

4.74 

4.81 

4.97 

4.63 

3.85 

2^2 

5.07 

5.18 

5.38 

5.03 

4.62 

B2 

5.66 

5.77 

5.89 

5.47 

5.12 

2^2 

6.26 

6.35 

6.51 

6.18 

5.52 

Bi 

6.53 

6.66 

6.83 

6.50 

5.90 

Ai 

7.02 

7.17 

7.37 

6.96 

6.70 

s-tetrazine 

Bzu 

2.97 

3.08 

3.07 

2.67 

2.25 

-^U 

4.16 

4.30 

4.32 

3.98 

3.40 

-^U 

6.00 

6.13 

6.11 

5.72 

5.00 

Bsu 

7.11 

7.25 

7.27 

6.95 

6.34 

MAD 

0.60 

0.72 

0.78 

0.49 


RMS 

0.64 

0.76 

0.82 

0.53 


MD 

1.04 

1.25 

1.24 

0.97 



TABLE III. Excitation energies for a number of single ring 
azabenzenes computed with the d-aug-cc-pVDZ basis using 
MP2/6-311-I-G** geometries. Units are in electron volts (eV), 
experimental references can be found in Ref. 0. 


tation energies for both EOM-LCCD and EOM-LCCSD, 
with average deviations relative to EOM-CCSD of 0.25 
and 0.35 eV respectively (see the correlation figures [2l|3] 
and Table iryi) . The comparison of the linear EOM meth¬ 
ods against EOM-CCSDT-3 is less satisfactory with a 
MAD of 0.45 and 0.56 eV for EOM-LCCD and EOM- 
LCCSD respectively. 

However, these methods are not without merit should 
the ground state wavefunction need to be computed with 
LCCD or LCCSD. This is an important concern for 
weakly interacting systems or transition states where an 
MBPT(2) wavefunction is a poor approximation. With 
an acceptable standard deviation of 0.11 and 0.10 eV 
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MBPT(2) 

LCCD 

LCCSD 

relative to EOM-CCSD 



MAD 

0.13 

0.25 

0.35 

RMS 

0.15 

0.28 

0.36 

MD 

0.37 

0.55 

0.80 

relative to EOM-CCSDT-3 



MAD 

0.31 

0.45 

0.56 

RMS 

0.38 

0.53 

0.60 

MD 

1.19 

1.49 

1.37 

TABLE IV. 

Excitation energy 

comparison 

statistics for the 

Miilheim test set. Units are in electron volts (eV), EOM- 

CCSD and EOM-CCSDT-3 reference values 

were taken from 


Sous et al^ 


compared to EOM-CCSD, the excitation energies and 
relative state orderings obtained with the linear EOM 
methods are reliably consistent compared to canonical 
EOM-CCSD. This is further illustrated in Tables HI HIl 
and m where the systematic and consistent overestima¬ 
tion of the excitation energies is balanced by a nearly 
exact agreement in state ordering and relative state en¬ 
ergy differences (a MAD of 0.04, 0.04 and 0.05 eV for 
EOM-MBPT(2), LCCD and LCCSD respectively). For 
approximate methods, getting these relative properties 
correct is just as useful as quantitatively accurate exci¬ 
tation energies. 


V. CONCLUSIONS 

We have expanded the coupled-cluster similarity trans¬ 
formed Hamiltonian using general coupled-cluster per¬ 
turbation theory to obtain the arbitrary order coupled- 
cluster perturbation theory effective Hamiltonian given 
by Eq. [T71 The inclusive (the sum of all orders 0 to 
n) form of the Hamiltonian is inserted into the stan¬ 
dard electronic-excitation equation-of-motion theory to 
give the completely general EOM-CCPT (Eq. 1^51) . The 
result is a way of using the standard equation-of-motion 
theory to directly compute excitation energies that are 
consistent with a given CCPT ground state wavefunc- 
tion. 

By choosing the generalized many-body perturbation 
theory partitioning of the Hamiltonian, we re-derive the 
well known EOM-MBPT(2)2a (EOM-CCSD(2)21) equa¬ 
tions. An alternative source of amplitudes is obtained by 
short-circuiting the inclusive {(l)g\H^^^\(j)o) set of equa¬ 
tions. The resulting infinite-order amplitudes correspond 
to the linear CCSD (and CCD in the case of single excita¬ 
tions being neglected) expansion, and their use in the in¬ 
clusive similarity transformed Hamiltonian in place of the 
standard MBPT(2) amplitudes provides an equation-of- 
motion method perturbatively consistent with a LCCD 
and LCCSD reference wavefunction. 

These approximate equation-of-motion methods are 


benchmarked by employing the Miilheinii^^— and 
Yale^i^ small organic molecule data sets. We use the 
diverse Miilheim test set to benchmark our new methods 
against the canonical EOM-CCSD and EOM-CCSDT-3 
methods, while the Yale test set is employed in bench¬ 
marking against experimental spectra. Our methods are 
found to consistently over estimate excitation energies, 
relative to EOM-CCSD and EOM-CCSDT-3 by - 0.25 
and 0.5 eV respectively. The precise statistics and 
correlation plots can be found in Table |IV] and Fig¬ 
ures [5] and O This systematic overestimation compared 
to the complete EOM-CCSD suggests that the similar¬ 
ity transformed Hamiltonian (and thus the excited state 
spectrum) based on the linear CC wavefunction is over¬ 
approximated compared to the ground-state description 
leading to an increased separation. The very good ~ 0.1 
eV standard deviation compared with EOM-CCSD sug¬ 
gests that the predicted relative spectra will be much 
more accurate than the absolute excitation energies. This 
is supported by studying the benchmark values from the 
Yale set given in Tables HI HIl and Hill where the presented 
approximate EOM methods obtain relative state order¬ 
ings and energies to within ^ 0.04 eV of canonical EOM- 
CCSD. 
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